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Abstract. We prove the existence of a resonance free region in scattering by a strictly 
convex obstacle O with the Robin boundary condition d v u + ju\go — 0. More precisely, 
we show that the scattering resonances lie below a cubic curve Im£ = — S\(\z + C. The 
constant S is the same as in the case of the Neumann boundary condition 7 = 0. This 
generalizes earlier results on cubic poles free regions [3] , [8] , [22] obtained for the Dirichlet 
boundary condition. 

I. Introduction and Statement of Results 

Let O be a bounded open set in W 1 with a smooth boundary. For 7, a real-valued smooth 
function on dO, we define P^ to be the Laplacian on the exterior domain — Ak«\o realized 
with the Robin boundary condition d u u + ^u\oo = 0. The operator P^ is a self-adjoint 
operator with continuous spectrum [0, 00). The resolvent and R^\Q : = (P^ — £ 2 )~ 1 : 
L 2 (R n \ O) L 2 (R n \ O) is holomorphic in Im( > 0. The Green function, GW(C, x, y) is 
defined as the kernel of this resolvent: 

u(x)= [ G<?>\t,x,y)f(y)dy, feC™(R n \0), 

JR n \0 

(-A - ( 2 )u{x) = f(x), i6i" \ O, d v u{x) + 7(x)w(x) = 0, 16 dO. 

The Green function has a meromorphic extension in ( to the whole complex plane if n is 
odd, and the logarithmic covering space of C \ {0} if n is even. 

The scattering poles or resonances of P^ are defined as the poles of this meromorphic 
continuation and they enjoy interesting interpretations and properties. The real part of a 
resonance corresponds to a frequency of a resonant wave, and the imaginary part, to the 
decay rate of the wave. Consequently understanding of the separation of resonances from 
the real axis is related to the decay properties and long time behavior of waves. 

Resonance free regions near the real axis have been extensively studied since the work 
of Lax-Phillips [10] and Vainberg [25] where the presence of such regions was linked to 
propagation of singularities for the wave equation, and hence to the geometry of the ob- 
stacle O. Thanks to the work of Melrose, Ivrii, Sjostrand and Taylor on propagation of 
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Figure 1. The Resonance- Free Region 

singularities for boundary value problems, if O is nontrapping, that is, all reflecting rays 
escape to infinity, there are no resonances in the region lm( > — Mlog |C| + Cm for any 
M - see [9, Chapter 24], [13], [2 4] and references given there. When the boundary dO is 
real analytic, and the obstacle is nontrapping, the work of Lebeau [11] on propagation of 
Gevrey 3 singularities implies that the resonance free region is cubic, that is, there are 
no resonances in lm( > — C 1 C 1 1,/3 + C\ for some constants C ,Ci - see Popov [15] and 
Bardos-Lebeau-Rauch [3] . The above results have been typically stated in the case of the 
Dirichlet boundary condition, as they depend only on propagation of singularities, they 
hold for ?W, with 7 analytic, for the cubic poles free region. 

The case of stictly convex obstacles has been studied, first when O is the sphere, since 
the work of Watson [26] on electromagnetic scattering by the earth almost a hundred years 
ago. In that case, for the Dirichlet or Neumann problems scattering poles are given by 
zeros of Hankel function or their derivatives, respectively - see Stefanov [18] for a modern 
account and references. Since a convex obstacle is nontrapping the existence of a logarithmic 
resonance free region follows from the general results [13]. It was first established, in the 
star shaped case, by Morawetz-Ralston-Strauss [14]. 

A remarkable discovery was made by Harge-Lebeau [8] who showed that for a smooth 
boundary and the Dirichlet boundary condition, the resonance free region is cubic, that is, 
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the same result is valid as for analytic non-trapping obstacles. This result is sharp as was 
shown already in [3] where the analysis of Gevrey-3 singularities of the wave trace gave a 
string of of resonances near a cubic curve - see also [17] and [20, Example 4]. The argument 
of Harge-Lebeau was based on the complex scaling method of Aguilar-Combes [1] and 
Balslev-Combes [2] developed for boundary value problems by Sjostrand-Zworski [19], [21]. 
In [22] Sjostrand-Zworski gave a more direct proof of the cubic resonance free region for 
smooth strictly convex obstacles and established polynomial bounds on the number of 
resonances in neighborhood of the real axis. In [23] they proved an asymptotic law for the 
number of resonances in the cubic strips for obstacles whose boundary satisfies a pinched 
curvature assumption. 

In this paper, we prove that the results of [8] and [22] on the cubic resonance free region 
are valid for an arbitrary Robin boundary condition. The main result which follows from 
Theorem 2 below is 

Theorem 1. Suppose that O is strictly convex with a smooth boundary. Then there are no 
resonances of in the cubic region 



with Ki{y) the principal curvatures of dO at y, —([ the first zero of the derivative of the 
Airy function. 

Remark 1. The constant S here is optimal if the obstacle O is a ball and 7 is a constant, 
since we can get explicit expressions for the resonances by Hankel functions - see [18]. 
When 7 is not a constant function Theorem 1 is new even in the case of the sphere. 

In the case of Dirichlet or Neumann boundary condition, the optimal value for the con- 
stant S is given in [3] and [17] when DO is analytic and its geodesies satisfy certain condi- 
tions. The optimal constant S in other cases is still unknown. 

The basic strategy is similar to that in [22] but with some challenges provided by the more 
general boundary condition. First, we write the operator in normal geodesic coordinates 
with respect to the boundary and apply the method of complex scaling all the way to 
the boundary. For this part, we shall follow the approach in [19] and [21] and reduce 
the problem of resonances free region to the problem of localizing the spectrum of a non- 
selfadjoint operator. Then we shall use the global FBI transform on the boundary [27] to 
change the complex-scaled Laplace operator to an Airy-type ordinary differential operator, 



ImC>-S|C|3+L7, 
for some C > 0, depending on 7 and O . The constant S is given by 



(1.1) 
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at least near the boundary. In Section 3, we study the model Airy- type ordinary differential 
operator with general boundary condition. In Section 4, we will go back to the scaled 
operator and prove a certain lower bound for it - see Theorem 2. Finally, in Section 5, 
we use that lower bound to show the existence of a resonance free region stated in the 
main theorem. In the Appendix we review the complex scaling up the boundary [21] as 
additional care is needed when dealing with general boundary conditions. 

Acknowledgement. I would like to thank Maciej Zworski for encouragement and advice 
during the preparation of this paper. Partial support by the National Science Foundation 
grant DMS-1201417 is also gratefully acknowledged. 

2. Preliminaries 

2.1. The complex scaling method. In this section, we reduce the problem of resonances 
to the spectrum of a non-selfadjoint operator by the complex scaling method. This method 
was first introduced by Aguilar-Combes [1] and Balslev-Combes [2] in studying the contin- 
uous spectrum of Schrodinger operators and later proved to be a strong tool in the study 
of resonances. We shall follow the methods of Sjostrand and Zworski and apply it to the 
case of Robin boundary condition. For more details, see [19], [21] and [22]. 

Let O be a strictly convex, bounded open set in W 1 with smooth boundary, then d(x) = 
dist(x, O) is a smooth convex function in Moreover, d"(x) ^ and dimker d"(x) = 1, 

generated by x — y(x) where y(x) is the closest point to x on dO, so that d(x) = \x — y(x)\. 
At y(x), the exterior unit normal vector of dO is u(y(x)) = Vd(y(x)) = uZy^i ■ 

Near every xo G dO, we shall choose the normal geodesic coordinates (x f , x n ), x' G U C 
M n_1 , x n ^ for IR n \C where (U, x' = (x\, . . . , £ n -i)) is a local coordinates for dO centered 
at xo, x n = d(x). This normal geodesic coordinates is given by x = s(x') + x n u(s(x')) where 
s : U — > dO is the coordinate map for dO. In the global version, the normal geodesic 
coordinates identify W 1 \ O with dO x [0, oo) by 

R n \0 3 x ^ (y, x n ) G dO x [0, oo) 

where y = s(x J ) G dO. 

Under this coordinates, for x n small, we have the following expression for —h 2 A: 

-h 2 A = {hD Xn f + R(x', hD x ,) - 2x n Q(x', hD x ,) + 0{x 2 n {hD x ,f) + 0(h)hD x + 0{h 2 ) 

where we introduce the semiclassical parameter h > which we will later let tend to zero, 
R and Q are two quadratic forms which are dual to the first and second fundamental form, 
respectively. Thus the principal curvatures of dO at y — s(x') are the eigenvalues of Q 
with respect to R. 
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Now following [21] and [22] we introduce the family of complex scaling contours (0 < 
0<O O ) 

T e = { z = x + idf(x) : x e R n \ O} c R n \ O + iW 1 

where / : W 1 \ O — > K is a smooth function such that for x near O, f(x) = ^d(x) 2 , so 
f'(x) = d(x)d'(x). 

If we parametrize Tgbj x E M. n \0, then we can compute the principal symbol of — A|r e : 
pe(x, = ((1 + i9f"(x))-% (1 + ief'ix))-^} =a e - ib e 

where 

a e = ((1 - (9f"(x))%£), b e = 29(f"(x)H). | = (1 + (9f"(x)) 2 )-^. 

In normal geodesic coordinates, we can write the contours Tg as the image of 

U x [0,oo) 3 (x',x n ) ^ s(x') + (1 + iB)x n u{s{x')) G C n , 

locally for x n small. In the global version, by rescaling t = |(1 + z#)|x n , the contours are 
the image of 

dO x [0, oo) 3 (y, t) y + G C\ 

where : [0, oo) — > C is a smooth injective map such that \g' e \ = l,g(0) = 0,g{t) = t ,\+f 6 , 
for t near 0; g(t) = outside a small neighborhood of and 

1 

arg(l + i<p(6)) ^ arg#(t) < arg(l + i6), - arg(l +i<p(6)) ^ arg^'(t) < (l + i9). 

We shall choose ip(8) small enough such that Tg satisfy the conditions given in [21]. In 
particular, we shall work on the contour r = Tg for 9 = 9i such that g = ge 1 (t) equals to 
te^ for t near (later on we shall use t G [0, L~ l ] for L large enough). On this contour, 

-/i 2 A| r = -^—(hD t ) 2 + R(y, hD y ) - 2g(t)Q(y, hD y ) 

+ 0(t 2 (hD y ) 2 ) + 0(h)hDy tt + 0(h 2 ), 

which is elliptic in both semiclassical sense and the usual sense. 
For t small such that g(t) = te^t , 

-/i 2 A| r = e-'^{{hD t f + 2tQ{y, hD y )) + R(y, hD y ) 
+ 0(t 2 (hD y ) 2 ) + 0(h)hD yjt + 0(h 2 ). 

Let p be the principal symbol of — h 2 A\ r . We notice that from [21], if g = g 9l , = ip{6\) 
are chosen suitably, then p lies in the lower half plane and for every 5 > 0, there exists 
e > such that 

t ^ 5 =>• e ^ — argp ^ n — e. 
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Now we consider the boundary condition. In the normal geodesic coordinate (y, t) G 
dO x [0, oo), the Robin boundary condition becomes 

d t u + 7u|t =0 = 

Therefore in the scaled operator, we shall choose the boundary condition 

e~^d t u + 7«|i =0 = 

or 

d v u + ku\gr = 
where k : dT — > C is a smooth function. 

We shall define the scaled operator P = Pp } = -A| r : D(T) ->■ L 2 (r) where D(T) = 
{u G /f 2 (r)|9yM + fc-u|gr = 0} and regard P as an operator on IR n \C by the parametrization 
of T given above. According to [19], P has discrete spectrum in the sector e~ l (°' 2tp >(0, +oo). 
Moreover, we have 

Proposition 2.1. The resonances of in —ip < arg£ < is the same as the square 
root of the eigenvalues of P in —2ip < arg£ < 0. 

The proof of the theorem is based on the following deformation results. First, we have 
the lemma about non-characteristic deformations which is proved in [19]. 

Lemma 2.2. Let u C lR n be an open set, h : [0, 1] x oj 3 (t,y) h- > h(t,y) G C n be a smooth 
proper map such that 

(1) det(d y h(t,y))^0for all (t,y); 

(2) hif,-) is injective; 

(3) hit, y) = h(0, y) for y G u\K where K is a compact subset of u. 

We write T t = h({t} x u). Let P(x, D x ) be a partial differential operator with holomorphic 
coefficients defined in a neighborhood of h([0, 1] x oo) such that P\r t is elliptic for ^ 
t ^ 1. If Uq G $'(Tq) and Pt Uq extends to a holomorphic function in a neighborhood of 
h([0, 1] x oo), then w extends to a holomorphic function in a neighborhood of h([0, 1] x oo). 

Next we need to the following lemma providing that we can apply the deformation all 
the way to the boundary and it will satisfy the desired boundary properties. 

Lemma 2.3. Let u G C°°(R n \ O) satisfy (-A - X 2 ) k °u = 0,d a u\ dn = u a G C°°(dtt) in 
a neighborhood of xq. Then there exists a complex neighborhood W of Xq such that (1) u 
extends holomorphically to a function U in a complex open neighborhood of W fl U|e|^e H?/ 

(2) uq = U\r g is smooth up to dTg = dO; 

(3) (-A|r 9 - X 2 ) ko u e = 0,d a u e \ dTe = u a in T e n W. 
Moreover, we may replace W 1 \0 by any fixed T v with \rj\ < 9q. 
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This lemma was proved in [21]. We recall the proof with more details in the appendix. 
Combining these deformation results with the method in [19], we have the proposition. 

2.2. Semiclassical Sobolev spaces. In this section, we review some basic facts about 
the semiclassical Sobolev spaces, especially the estimate of the trace operator which we 
shall use later. 

Let H^W 1 ) C iS'(IR n ) be the semiclassical Sobolev space of order s with the norm 

IMItf£(R») = \\(hD) s u\\ L 2 {Rn) , (hD) = (1 + (hD) 2 ) 1 ?. 

Let X be a compact smooth manifold, we can choose a finite cover X\ , . . . , Xp of X 
where X±, . . . , X p are coordinate charts with local coordinates xi,...,x n . Then there exists 
a partition of unity Xj C£°(Xj), Y7j=i Xj = 1- We define the semiclassical Sobolev space 
HfJ^X) to be the space of all u G D'(X) such that 

v 

IMIfljjW = Yl \\Xj( hD ) s Xju\\l2 (Xj) < oo. 
i=i 

For different choice of the coordinate charts and partition of unity, the norms are equiv- 
alent uniformly for h > 0. Also, another equivalent norm can be given by ||u||h s (X) — 
II (/ - h 2 A) 2 u\\l 2 {x)i where A is the Laplacian operator with respect to some Riemann- 
ian metric. From this norm, we see that H^(X) is a Hilbert space with inner product 
(u,v) H s h(x) = ((I - h 2 A)lu, (7 - h 2 A)h) L 2 {x) . 

Now consider the trace operator Tr : C°°(lR n \ O) — > C°°(dO),u i— >■ u\qo- In the normal 
geodesic coordinates given in the previous section, it is equivalent to the operator Tr : 
C°°(X x [0,oo)) -»■ C°°{X),Tru{y) = u{0,y). 

Proposition 2.4. For u G Cfi°(X x [0, oo)), we have || Tr < ^ -1 IMIff 2 (xx[o oo))- 

Proof. Since u G C^°(X x [0, oo)), we know there exists L > such that u is supported in 
X x [0,L]. Therefore 

\\^AHl<iX) = - h ~ 1 / hD t \\u(;t)\\ 2 Hl{x) dt^2h- 1 / \(hD t u(;t),u(;t)) H i ix) \dt 
</ J 

/*oo 

< h- 1 / [\\hD t u{.,t)\\V h{ x) + l|t*0,*)ll^ W ]* < C^ 1 |KII^(xx[o I oo))- 

** 

□ 

Remark 2. In fact, we can even prove that Tr = O s [h 2 ) : H^(ft) — > H h 2 (dQ) when 
s > | ; Q CC lR n an open set with smooth boundary. 
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2.3. The FBI transform. This section is devoted to a brief introduction of Fourier-Bros- 
Iagolnitzer (FBI) transform on IR n or a smooth compact manifold. The FBI transform was 
first introduced by Bros and Iagolnitzer to study the analytic singularity of a distribution 
and later on became an extremely useful tool in microlocal analysis. We refer to the books 
by Delort [5], Folland [7] and Sjostrand [16] for more details. An account of the semiclassical 
theory needed here can be found in Martinez [12] and Zworski [28, Chapter 13]. In our 
argument an important part is played by the proof of the sharp Garding inequality given 
by Cordoba and Fefferman [4]. Wunsch and Zworski [27] (see also [17] for the analytic case) 
adapted the FBI transform to compact Riemannian manifolds. 

We shall first review the basic facts about FBI transform on IR n , then follow [27] to 
describe the FBI transform on a compact manifold which for us will eventually be DO. 

The FBI transform of a function u E J^ / (R n ) is given by 

T h u(x,£) = 2-?(nh)- § r [ e-^ x - y)2+ ^ x - yH u(y)dy. 

This is a special case of the wave packet transform 

T h u{x,£) = (u,(p( X £,h))y,y, 

where <P( X £,h) is a "wave packet" concentrated at (x,£) E T*IR n . If we take <f( x ,£,h) to be the 
coherent state, then we get the FBI transform. The FBI transform is also closely related 
to the Bargmann transform in several complex variables: 

f h u{z) = [ e-^ {z - y)2 u(y)dy,z E C n . 

If we identify z = x - i£ E C n with (x, f ) E T*R n , then 

, .v n , . , 3n 1 /-2 ~ t ■. 

T h u(x,0 = 2-2{irh)~e-^ T h u{z). 

Let us define L|(C n ) to be the L 2 -space with the weight e'^^^m^dz) where m(dz) is 
the Lebesgue measure on C n , then Th : L 2 (M. n ) — > L|(C n ) is an isometry with image 
#$(C n ) = {f E L|(C n ) : / is holomorphic in C n }. 

Back to the FBI transform T h : L 2 (R n ) -> L 2 (R 2n ), the adjoint T* h : L 2 (R 2n ) -> L 2 (R n ) 
is given by 

T*v(y) = 2-f (nh)~^ [ e-^ x - y)2 -^ x - yH v(x^)dxdi 

Ji 2 " 

From the properties of the Bargmann transform, we know that Th is an isometry with 
image L 2 (R 2n ) n e~^A(C™_^). Thus T*T h u = u for every u E L 2 (R n ) and T h T* is the 
orthogonal projection in L 2 (IR 2n ) onto the image L 2 (R 2n ) D e~ ^ A(C^_^) . 
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Now let (X, g) be a compact Riemannian manifold. We write y to be a point on X, dy 
to be the volume form on X; (x,£) a point on T*X (where G T*X) and dxd^ to 

be the canonical volume form on T*X. 

An admissible phase function <p(x,£,y) is a smooth function on T*X x X satisfying the 
following conditions: 

(1) if is an elliptic polyhomogeneous symbol of order one in £; 

(2) kmp ^ 0; 

(3) d y (f\ A = -£dy; 

(4) 4lm^| A ~<£); 

(5) <^|a = 0. 

Therefore near the diagonal A, <p = £ ■ (sc — y) + (Q(x, £, y)(x — y),(x — y)) where Q is a 
symmetric matrix- valued symbol of degree 1 in £ with IulQ|a ~ (0-^- 

The symbol class h m S^ hg {T* X x X) is defined to be the collection of all smooth functions 

a = a(x, £, y; h) ~ h m (a k (x, f , + ha k ^(x, £,y) H ) 

where a,j(x,£,y) are polyhomogeneous symbols of degree j in £ and the asymptotic expan- 
sion means that 

\a - h m (a k + ■■■ + tf'a fc _i)| ^ C 3 h m ^ +1 \^-\ |£| > 1. 

The symbol class h m Sp h (T* X) is defined similarly without the y-components. 

A symbol a G h m Sp hg is called elliptic if the principal part \a k \ ~ (£) uniformly with 
respect to other variables. The quantization of a is defined as the operator 

Op(a)u(y) = y e^^^a^^y; /i)u(x) X (y,x)dxde, 

where exp is the exponential map with respect to g> on X. We shall write h m ^ h (X) to be 
the algebra of all operators Op(a) +fi,a6 /i m S fc , = 0{h°°) : C-°°(X) -> C°°(X). 

Then following [27], an FBI transform on X is an operator T h : C°°(X) ->■ C°°(T*X) 
given by 

T h u(x,0= [ e^^a(x,C,y,h)x{x,^y)u(y)dy. (2.1) 

Here (p(x,£,y) is an admissible phase function; a(x,£,y; h) G h^^S^ is an elliptic polyho- 
mogeneous symbol; ^(a;, 2/) is a cut-off function to a small neighborhood of the diagonal 
A = {(x,£,y) G T*X x X : x = y} such that Im^ ^ —C~ 1 d(x,y) 2 on the support of \. 
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The following properties of the FBI transform were proved in [27]: 

(1) T h : L 2 (X) ->■ L 2 (T*X) is bounded for h < h ; 

(2) We can choose a suitable phase if and elliptic symbol a such that 

\\T h u\\ LHT , x) = (1 + 0(/0)NU*(X), (2.2) 

i.e. T/j is an isometry modulo h°°. 

From now on, we shall always use such kind of FBI transforms. Furthermore, we know 
from [27]: 

Lemma 2.5. Let P = Op(p) G h k ^ m (X), then T*pT h - P G h k+1 ^ m -\ 

We can apply this to prove the following 
Proposition 2.6. (1) For any u G C°°(X), 

\\(0 T hU\\L2(T*X) ^ C\\u\\ H i {x) . (2.3) 
(2) If A(x, hD x ) is a second-order differential operator on X, then for any u G C°°(X), 
\\A(x,0Tu\\ 

L 2 (T*X) = \\A(x,hD x )u\\ 2 L H x) + 0(h)\H\ 2 Hl (2.4) 

Proof. (1) 

\\(0T h u\\ 2 = ((0T h u,(0T h u) = (T*(0 2 T h u,u) 

= ((I - A)u,u) + (Ru,u) = + (Ru,u), 

where R G h®\ So (Ru,u) = 0(h)\\u\\ 2 , . 

(2) Notice that by symbol calculus 

(AA)(x, hD) = A(x, hD)*A(x, hD) mod h^ 3 

we have the following 

\\A(x,0Tu\\ 2 LHT , x) = (A(x,OTu,A(x,£)Tu) = (T*\A(x,0\ 2 Tu,u) 

= (A(x, hD)*A(x, hD)u, u) + (Ru, u) = \\A(x, hD)u\\ 2 + (Ru, u), 

where R G h^ 3 . So (Ru,u) = 0(h)\\u\\ 2 3 • □ 

Remark 3. We also notice that all the discussion above work for functions with value in 
a Hilbert space J^. In our case, we shall choose X = dO and <ffl = L 2 ([0, oo)). 



resonance-free region in scattering by a strictly convex obstacle 11 
3. Estimates of the Airy- type operator 

In this section, we shall give the lower bounds for the ordinary differential operator 

P = e -2«/3(( /ijDt )2 + t j + (h)hD t + 0(h + hht + t 2 ) (3.1) 

define on [0, oo) with general conditions at t — 0. 

3.1. The Dirichlet and Neumann realization. Let Ai be the Airy function defined by 

T 3, 



Ai(s) = ^z f e< a3 M +ias d(j. 

° 'lmo-=5>0 



2tt 



Then we can give all the eigenfunctions and eigenvalues for the Dirichlet and Neumann 
realization of the Airy operator D 2 + s on [0, oo): 

(D 2 S + s)Ai(s - Q = QAi(s - Q, Ai(-Q) = 0; 

(D 2 + s)At(s - Q = C 3 Ms - Q, Ai'(-Q = 0, 

where < Ci < C2 < ■ • ■ , < ([ < ('2 < ■ ■ ■ and Ci ~ 2.338, ([ « 1-019. The spectral 
theorem gives the following estimates: 

Let v E Q°[0, 00), if u(0) = 0, then 

((D 2 8 +s)v,v) >Ci\\v\\ 2 ; (3.2) 

if D s v(0) = 0, then 

{(D 2 s +s)v : v)>av\\ 2 - (3-3) 
Now we consider the semiclassical version of the Airy operator, (hD t ) 2 + 1. By changing 

2 

the variables t = h~ss, we have 

(hD t ) 2 + t = hl(D 2 + s). 

12 12 

For -u = u(t), we define u(s) = h3u(h3s), then = h~3v(h~3t), \\u\\ L 2 = \\v\\ L 2 and 

((hD t ) 2 + t)u(t) = hi-(D 2 + s)v(s). 

Therefore 

(((/lA) 2 + *K«)xf = h 2 H(D 2 + S )v) L , 
Applying the estimates (3.2) and (3.3), we have for u 6 C£°[0, 00), if u(0) = 0, then 

(((hD t ) 2 + t)u,u) > Ci^NI* (3.4) 

if D t u(0) = 0, then 

{((hD t ) 2 + t)u,u)>([h*\\u\\ 2 (3.5) 
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We also have the following useful identity: for u G Cq°([0, oo)),w(0) = or D t u(0) = 0, 
{{{hD t f + t)u,u) = ((hD t ) 2 u,u) + (tu,u) = \\hD t u\\ 2 + ||t^|| 2 , 
and consequently 

(((hD t ) 2 + t)u,u) > \\hD t u\\ 2 ; (3.6) 

(((hD t ) 2 + t)u,u)^\\t-*u\\ 2 . (3.7) 
Now we could estimate \\((hD t ) 2 + t)u\\ by the Cauchy-Schwartz inequality 

\\((hD t ) 2 + t)u\\\\u\\ > {{{hD t ) 2 +t)u,u). 

If u(0) = 0, then by (3.4) 

\\((hD t ) 2 + t)u\\ >(ihs\\u\\, (3.8) 

and by (3.6) 

\\{{hD t f + t)u\\ 2 > C,hl {{{hD t f + t)u, u) > C,hl \\hD t u\\\ 

thus 

\\((hD t ) 2 + t)u\\ > y/&h*\\hD t u\\ (3.9) 
Similarly, if D t u(0) = 0, by (3.5) and (3.6), we have 

\\((hD t ) 2 + t)u\\ ^ Ci^lMI; (3-10) 

\\{(hD t f+t)u\\ > ^C[h^\\hD t u\\. (3.11) 

Another way to estimate \\{{hD t ) 2 + t)u\\ is to use the following identity: If u(0) = or 
D t u(0) = 0, since (u, hD t u) is real 

\\((hD t ) 2 + t)u\\ 2 = \\{hD t fu\\ 2 + INI 2 + 2 Re(tu, (hD t ) 2 u) 

= \\(hD t ) 2 u\\ 2 + \\tu\\ 2 + 2Re{hD t (tu), hD t u) 

= \\{hD t ) 2 u\\ 2 + \\tu\\ 2 + 2 Re(thD t u, hD t u) + hRe-(u, hD t u) 

= n(/iA) 2 «ii 2 +iNr+2iit2/ l D t tt|i 2 . 

This gives us the following estimates 

IK^Af + ^ir^ll^A^f + II^II^IK/iA) 2 ^! 2 . (3.12) 
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3.2. General condition. Now we remove the Dirichlet or Neumann condition at t = 
and try to get a lower bound of (((hD t ) 2 + t)u,u). In this case, (hD t ) 2 + t is no longer a 
self-adjoint operator, but the semiclassical setting allows us to view it as a perturbation of 
the Neumann realization. We shall start with the following basic estimate: 

\\hD t u\\ 2 = (hD t u,hD t u) = ((hD t ) 2 u, u) - ih 2 D t u(0)u{0) 
sC (((hD t ) 2 + t)u, u) - ih 2 D t u(0)u(0). 
Since the right hand side is real, we have 

\\hD t u\\ 2 sC Re(((hD t ) 2 + t)u,u) - Re(ih 2 D t u(0)u(0)) 
sC Re(((hD t ) 2 +t)u,u) + h 2 \D t u(0)\\u(0)\ 
^ Re(((hD t ) 2 + t)u,u) + O(h 2 )\D t u(0)\ 2 + O(h 2 )\u(0)\ 2 , 

or 

Re(((hD t ) 2 + t)u,u) > \\hD t u\\ 2 - O(h 2 )\D t u(0)\ 2 - O(h 2 )\u(0)\ 2 . (3.13) 

which is the analogue of (3.6) for general u. Next we try to get an analogue of (3.4) and 
(3.5). 

Lemma 3.1. Suppose u G Cq°([0, oo)) ; then we have the following estimate: 

Re(((hD t ) 2 +t)u,u) ^ C[hHl - 0(hl))\\u\\ 2 - O(h 2 )\D t u(0)\ 2 , (3.14) 

| lm(((hD t ) 2 + t)u, u)\ ^ 0(0) Re(((hD t ) 2 + t)u,u) + O(h 2 )\D t u(0)\ 2 . (3.15) 

2 12 

Proof. Write D t u(0) = a for simplicity. First, we use the scaling t — h^s, v(s) = hzu(hzs) 
as before. We have 

(((hD t ) 2 + t)u,u) L 2 = hl((D 2 s + s)v,v) L 2, 

and more importantly, D s v(0) = hD t u(0) = ha. Now let v(s) = w(s) + has\(s) where 
X € Co°([0, oo)) is a fixed function such that x = 1 near 0. Then we get the decomposition 

((D 2 + s)v,v) = 

((D 2 + s)w, w) + ha((D 2 s + s)s X , w) + ha((D 2 + s)w, s X ) + h 2 \a\ 2 ((D 2 + s)s X , s X ). 

Since w satisfies the Neumann condition D s w(0) = 0, we know the first term is real. 
We can integrate by parts to rewrite the third term as —ha(D s w,D s (sx)) + ha(w, s 2 x) ■ 
Therefore for the real part, by the Cauchy-Schwarz inequality, 

Re((D 2 s + s)v,v) ^ ((D 2 + s)w,w) - 0(/i)|o|||D 8 u;|| - 0(h)\a\\\w\\ - 0(h 2 )\a\ 2 

^ ((D 2 + s)w,w) - 0(hl)\a\ 2 - O(0)(\\w\\ 2 + \\D s w\\ 2 ). 
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From 

\\D s w\\ 2 = (D s w,D s w) = (D 2 s w,w) ^ ((D 2 s +s)w,w) 
\\w\\ 2 KC 1 - 1 ((D 2 a + s)w 1 w) (by(3.3)), 

we have 

Re((D 2 s + s)v,v) > (1 - 0(hl))((D 2 + s)w,w) - 0(hl)\a\ 2 . (3.16) 

By (3.3) and 

Ml 2 > (IMI - \\has X \\) 2 = \\v\\ 2 - 0(h)\a\\\v\\ + 0(h 2 )\a\ 2 
^ (1 -0(hl))\\v\\ 2 -0(hl)\a\ 2 , 

we get 

Re<(£>2 + s)v,t;) ^ Ci(l - 0(/i*))H| 2 - 0{hl)\a\ 2 
>&(l-0{hi))\\v\\ 2 -0{hi)\a\ 2 . 

For the imaginary part, we have 

| Im((L> 2 + s)v,v)\ < 0(/it)|a| 2 + 0(hl)(\\w\\ 2 + ||D s w|| 2 ) 
< 0(hl)\a\ 2 + 0(hl)((D 2 s + s)w,w). 

Using (3.16) again, we have 

\lm({D 2 s +s)v,v)\ s= 0{hl)\a\ 2 + 0{hl)Re({D 2 s + s)v,v). 
Scaling back from v to u, we get the desired estimates (3.14) and (3.15). □ 

Finally, we need an analogue of (3.12). The argument in Section 3.1 shows that 

\\((hD t ) 2 + t)u\\ 2 = \\(hD t ) 2 u\\ 2 + \\tu\\ 2 + 2\\thD t u\\ 2 + hRe-(u, hD t u). 

% 

Although the last term is no longer zero, we can calculate it through integration by parts. 
Since 

(u,hD t u) = (hD t u,u) — hi\u(0)\ 2 , 

we have 

lm(u,hD t u) = -^|m(0)| 2 . 

Therefore 

2 

Re-(w, hD t u) = -%(0)| 2 , 

and we get 

\\((hD t ) 2 + t)u\\ 2 > \\(hD t ) 2 u\\ 2 - h 2 \u(0)\ 2 (3-17) 
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3.3. Restriction to a small interval. Now we restrict the support of u to a small fixed 
interval and get a better estimate. 

Lemma 3.2. If L > is sufficiently large, < h < h {L), then the following estimates 
holds uniformly for u G C^°{[0, L^ 1 ]): 

Re(((hD t ) 2 + t)u,u) ^ (C[hl - 0{hL))\\u\\ 2 - O{h 2 )\D t u{0)\ 2 + ^\\tu\\ 2 . (3.18) 
Proof. Choose Xo, Xi £ C°°(R) such that 

Xo = 1 on (-oo,/i2] 
Xi = 1 on [2hz, oo) 
<9 a x, = O Q (/rf),a = 0,l,2. 

Then we can deduce that 

Xo[xo, (hD t ) 2 ] + xi[xi, (hD t ) 2 ] = -[ Xo {hD t ) 2 { X o) + Xi(^A) 2 (xi)] = 0(h), 
from which we have 

({{hD t ) 2 + t)u, u) = ( X o{{hD t ) 2 + t)u, Xo u) + ( X i{{hD t ) 2 + t)u, X iu) 

= {{{hD t ) 2 + t)xoU, Xou) + {{{hD t ) 2 + t) Xl u, Xl u) 

- ((Xo[Xo, {hD t ) 2 } + Xi[Xi, {hD t ) 2 ])u, u). 

Since Xi M (0) = 0, also X $u and u have the same condition at t — 0, we can apply the 
estimates (3.7) and (3.14) to get 

Re({{hD t ) 2 + t)u,u) > Ci^ f (l - 0{hl))\\ X ,u\\ 2 + ||^xiw|| 2 - O{h 2 )\D t u{0)\ 2 - 0{h)\\u\\ 2 . 
From the construction of Xo, Xl> it is eas Y to see 

\\ u \\ 2 = \\Xou\\ 2 + WxMl 2 , \\tu\\ 2 = \\t Xo u\\ 2 + ||txiwf • 
Therefore we only need to prove for some cq = 0{L), we have 

c h\\xou\\ 2 > ^\\t X ou\\ 2 (3.19) 

and 

\\thiu\\ 2 > {(hi -c h)\\ Xl u\\ 2 + ^\\t Xl u\\ 2 . (3.20) 

Since xo is supported on [—00,2/12], we only need to choose Co = 2L to get (3.19). To 
prove (3.20), we need to show that for t G [h?, L^ 1 ], t ^ (hs — c oh + \t 2 or equivalently, 

L, ._,. . 1 



2 {t-L- 1 ) 2 ^— + 2hL-( 1 ht. 
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The left hand side achieves its maximum at t = hh ^ so we only need 

hk ^ (0 - c h + ^-h 

which can be achieved by choosing h < ho(L) small enough. □ 

3.4. Airy operator with lower order terms. Now we shall include the lower order 
terms and prove the main result of this section. 

Proposition 3.3. Suppose that a second order ordinary differential operator on [0, oo) 
satisfies that 

P = e -^^hD t f + t) + 0(h)hD t + 0(h + hh + t 2 ) 

Let Uq G C with argoj G ("I 71 ") I 71 ")- If L > is sufficiently large, h > sufficiently small 
depending on L and Uq, then for u G Cq°([0, L" 1 )), 

\\(P - co )u\\ 2 ^ (\e 2 ^u - C^ f i 2 - 0(hL))\\u\\ 2 2 
+ C^WihDtfuf - 0(h 2 )\D t u(Q)\ 2 - O(h 2 )\u(0)\ 2 . 

Proof. We begin with the following identity 

\\(P-lo )u\\ 2 = \\(e~^((hD t ) 2 + t)-u )u\\ 2 + \\[O(h)hD t + O(h + hh + t 2 )]u\\ 2 
+ 2 Re((e-^((/iA) 2 +t) - u )u, [0(h)hD t + 0(h + hh + t 2 )]u) 
> || (e^((/iA) 2 +t)- u )u\\ 2 - [0(h)(((hD t ) 2 + t)u, hD t u) 
+ 0(h)(u, hD t u) + (((hD t ) 2 + 0(/i + t 2 )u) + («, 0(ft + t 2 )u)}. 

The lower order terms are estimated as follows, 

0(h) (((hD t ) 2 + t)u, hD t u) ^ 0(h) \\((hD t ) 2 + t)u\\ 2 + 0(h) \\hD t u\\ 2 

^ 0(h) \\((hD t ) 2 + t)u\\ 2 + 0(h) Re(((hD t ) 2 + t)u, u) 
+ O(h 3 )\D t u(0)\ 2 + O(h 3 )\u(0)\ 2 (by (3.13)); 
0(h)(u,hD t u) <: 0(h)\\u\\ 2 + 0(h)\\hD t u\\ 2 

^ 0(h)\\u\\ 2 + 0(h) Re(((hD t ) 2 + t)u, u) 

+ O(h 3 )\D t u(0)\ 2 + O(h 3 )\u(0)\ 2 (by (3.13)); 
(((hD t ) 2 + t)u, 0(h + t 2 )u) ^ 0(h)(((hD t ) 2 + t)u,u) + 0(l)(((hD t ) 2 + t)u,t 2 u) 

^ 0(h)\\((hD t ) 2 + t)u\\ 2 + 0(h)\\u\\ 2 
+ l\\((hD t ) 2 + t)u\\ 2 + 0(l)\\t 2 u\\ 2 , 
(u,0(h + t 2 )u) 0(/i)||w|| 2 + 0(l)|H| 2 . 
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For the leading terms, we use the following identities 

\\(e-^((hD t f +t)- co )u\\ 2 = \\((hD t f + t)u\\ 2 + \u Q \ 2 \\uf 

- 2Re(e~^ ((hD t ) 2 + t)u,u u), 

and 

-2Re(e"^((/iA) 2 + t)u,u u) = 2 Re[e^u (((hD t ) 2 + t)u, u)} 

= Re(2e^a;o) • Re({{hD t ) 2 + t)u,u) 
- Im(2e^cJ ) ■ lm(((hD t ) 2 + t)u, u). 
By (3.15), the second term is bounded below by 

-2|w ||Im(((/iA) 2 +t)u,u)\ ^ -0(hl)Re(((hD t ) 2 + t)u,u) -O(h 2 )\D t u(0)\ 2 . 
Therefore 

- 2 Re(e"^ ((hD t f + t)u, u u} 

> (2cos(f -argw ) - O(h^))\uj \Re(((hD t ) 2 + t)u,u) - O(h 2 )\D t u(0)\ 2 . 
Now combining all the terms together, we get the following estimate 
\\(P - uj )u\\ 2 ^(2cos(f -argwo) - O(hl))\u \Re(((hD t ) 2 + t)u,u) 

+ (\u \ 2 - o{h))\\ u \\ 2 + (|- o(h))\\((hD t y + t)u\\ 2 

- 0(l)\\tu\\ 2 - 0(l)\\t 2 u\\ 2 - O(h 2 )\D t u(0)\ 2 - O(h 3 )\u(0)\ 2 . 

Since || — argw | < f , cos(| — argwo) > 0, when h is small enough, we can apply (3.18) 
to the first term and use 

\u \ 2 - 2Re(e-^cJ )Ci^ = \e 2m/3 co - Ci/^| 2 - 0(h) 

to get 

1 2 



||(P - u;oH| 2 > ((|e^w - T - C(/iL)) || W || 2 + (§ - 0(/i)) ||((/*A) 2 + t) 
+ (|w | cos(f - argwo) - 0(h))L - 0(1)) \\tu\\ 2 
- 0(l)||t 2 w|| 2 - O(h 2 )\D t u(0)\ 2 - 0(h 3 )\u'^ 2 



u 



Since u is supported in [0, L 1 ], we have ||t 2 w|| 2 ^ L 2 \\tu\\ 2 . So if L > L (u ) large 
enough, h < h (L,u ) small enough, we get 

|| (P - u )u\\ 2 > (\e 2m / 3 u - C(^| 2 - 0(hL))\\u\\ 2 + (| - 0(h)) \\((hD t ) 2 + t)u\\ 2 

+ C-'LWtuW 2 - O(h 2 )\D t u(0)\ 2 - O(h 3 )\u(0)\ 2 . 

Applying (3.17), we conclude the proof of (3.21). □ 
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Remark 4. If we replaced (hD t ) 2 + t by (hD t ) 2 + Qt, then the estimates (with ([0 replaced 
by ([Qshs ) remain uniform for Q G [C~ l ,C], \u \ G [C _1 ,C] and arg(w ) G [—^ + 5,^-— 5] 
(S,C>0) 

Remark 5. For the Dirichlet and Neumann realization, we can get the same estimate with- 
out the last two lower order terms — O(h 2 )\D t u(0)\ 2 — O(h 2 )\u(0)\ 2 based on the inequalities 
in Section 3.1. For Dirichlet realization, we can also improve C,[ to £i- 

4. Lower bounds on the scaled operator 

In section 2.1, we denned the scaled operator P = — A|r and find the explicit formula 
in normal coordinates with respect to the boundary: 3 x i— y (y, t) G X x (0, oo). If we 
freeze (y, 77) G T*X, then the (semiclassical) symbol of h 2 P is given by 

P(y, t; 77, hD y ) = e"T ((hD t ) 2 + 2tQ(y, 77)) + R(y, 77) + 0(t 2 + /i) (r/) 2 + 0(h)hD t . (4.1) 

In this section, we first estimate P(y, t; 77, hD t ) — ujq and then through the FBI transform 
introduced in Section 2.3 to get a lower bound on h 2 P — ujq. 



4.1. Estimate in the glancing region. When Rewo — R(y,rj) is small, the main part of 
P(y, t; 77, hD t ) — ^0 is given by the Airy-type operator e _ ~? ((hD t ) 2 + 2tQ(y, 77)). We shall 
apply our results in Section 3 to get the following estimate. Notice that such (y, 77) lies in 
a compact subset of T*X. 

Lemma 4.1. Let ujq G C with Reu; > 0,lmw = r > 0. Suppose \ Reu — R(x',^')\ < c 
where c is sufficiently small, L is large enough and < h < h (L). Then For any v G 
C^QO,^ 1 )), we have 

\\{P{y,t,r],hD t )-u )v\\ 2 > |r + 2S(Reu )h* - 0(h)\ 2 \\v\\ 2 

+C- 1 \\(hD t ) 2 v\\ 2 - O(h 2 )\D t v(0)\ 2 - O(h 2 )\v(0)\ 2 

where S is given by 



(4.2) 



2 

7T N 



S = 2-s cos ( - ) Ci ( . _ min ^ K^y) ) . (4.3) 



l\~ :( II) \ ' 

.6/ \yeX,j=l,...,n-l 

Kj(y) are the principal curvatures of X = DO at y and —([ is the first zero of the derivative 
of Airy function. 



Proof. Let c = tq tan |, then since | Re Uo—R(y, rf)\ < c, we have arg(co> — R(y, rf)) G [ 
It follows immediately from (3.21) by replacing uo with uq — R(y,rj) that 

\\(P(y,t,r),hD t )-u )vf > (\u - R{y,rj) - e~^ Ci(?Q(y^)) l ^\ 2 - 0(hL))\\v\\ 2 

+C- 1 \\{hD t ) 2 v\\ 2 - O(h 2 )\D t v(0)\ 2 - 0{h 2 )\v(^ 2 



7T 27T] 

3' 3 J 
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The uniformity of the constants follows from the ellipticity of Q and R. 
Now we need to find a uniform lower bound for 

\u - R(y,r)) - e"T C((2Q(y,r ? ))I/i§| 2 = | Reu - R(y,rj) + sin (|) c^Oi//. /,))"// : 



+ r + cos . . 

.0/ 



6. 

7T 



) Ci(2Q(2/,77))l^ 



(4.4) 



over (y,rj) such that |Rewo — R(y,r])\ < c. The minimum is obtained at R(y,rj) = 

2 2 

Keuo + O(h^) and the minimum of ([(2Q(y, r/))s under such constraint. Since the principal 
curvatures are the eigenvalues of the quadratic form Q(y, rj) with respect to the quadratic 
form R(y,rj), we have 

Q(Viv)> ( v m i n , K j(v) ) R {y,v)- 

\y£X,]=l,...,n-l I 



Thus 

(4.4) ^ \r + 2S{Reu )hl\ 2 + 0{hl) 
which completes the proof of (4.2). □ 

4.2. Estimate away from the glancing region. When |Reco>o — R(y,rj)\ is bounded 
from below, 2tQ(y, rj) is dominated by R(y, rj) — u for small t. In this case, we can give a 
better estimate for P from e~ir {hD t ) 2 + R(y,rj) — u . 

Lemma 4.2. Suppose that ujq G C with Rec^o > 0,Imuo = r > 0, \ Reuo — R(y,T])\ > c, 
then for L large enough, h sufficiently small, and v G C^°([0, L~ x )), 

\\(P(y,t,r), hD t ) - u )v\\ 2 > (r + C^fWvf + C- x {\\{hD t fvf + (v) 4 \\v\\ 2 ) 

-O(h 2 )(ri) 2 \v(0)\\D t v(0)\- 

Proof. Since 

[P{y, t, n, hD t ) - uj ]v = \e- 2 ^{hD t ) 2 + R(y, rj) - u ]v + [0(t 2 + t + h)(r]) 2 + 0(h)hD t ]v, 
we have 

\\(P(y,t,ri,hD t )-uo)v\\ 2 > (|| [e^(hD t ) 2 + R(y, rj) - u ]v\\ 

-\\[0(t 2 + t + h)(r,) 2 + 0(h)hD t }v\\) 2 . 

Now 

\\[e-^(hD t ) 2 + R(y,r))-cj }v\\ 2 = \\(hD t ) 2 v\\ 2 + \R(y, V ) - cu \ 2 \\v\\ 2 

2 ( 4 - 7 ) 

+ 2 Re[e-^{R(y, rj) - co )((hD t ) 2 v, v)\, 
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where 

Re{e-^{R{y,r]) - L0 ){{hD t ) 2 v,v)\ = Re[e-^{R{y,r]) - Q )] Re{{hD t ) 2 v,v) 

- lm[e'^ (R(y, rj) - Q )] lm{(hD t ) 2 v , v). 

Notice that 



((hD t ) 2 v,v) = \\hD t v\\ 2 + ih 2 D t v(0)v(0). 

Therefore 

Re((hD t ) 2 v,v) = \\hD t v\\ 2 + Re(ih 2 D t v(0)v(fij) ^ -h 2 \D t v(0)\\v(0)\; (4.9) 



\Im((hD t )\v)\ = \ Im(ih 2 D t v(0)v(0))\ < h 2 \D t v{0)\\v(0)\. (4.10) 
We can compute that 



2 1 s/3 

Re [e - ^ (R{y, rj) - cJ )] = ~-{R{y, rj) - Reu ) + —r Q . 



In the identity 



Re[e-^(i2(y,77) - u) )] Re((/iA) 2 ^) =- \{R{y,rj) - Reu )Re((hD t ) 2 v,v) 



+ ^-r Re((hD t ) 2 v,v), 



we apply (4.9) to the second term and 



\{R{y,rj) - Rec^o) Re((hD t ) 2 v,v)\ < ^(||(/iA) 2 ^|| 2 + K^M) - Re c^ ) 1 2 II ^ II 2 ) 



to the first term, and get 

Re[e-^(R(y,r)) - u )]Re{(hD t ) 2 v,v) > - h\\hD t v\\ 2 + \R(y, V ) - Reu \ 2 \\v\\ 2 ) 

(4.11) 



-/i 2 ^r |AT;(0)|K0)|. 
By (4.10), we also have 

-lm[e-^(R(y, V )-cu )]lm((hD t ) 2 v,v) > -h 2 \R(y,r]) - u \\D t v(0)\\v(0)\. (4.12) 
Combining (4.7), (4.8), (4.11), (4.12) together, we have 

\\[e-^{hD t f + R(y :V ) - u ]v\\ 2 > \^hD t fv\\ 2 + (r 2 + \\R{y^) - Reu, | 2 )|M| 2 

- h 2 (\R{y,rj) - o; 1 + V3r o )|A«(0)|K0)|. 

Now we estimate the remainder terms, since \R(y,rj) — Rec^ol > c and R(y,rj) is a 
quadratic form in 77, we have 

Cr 1 ^) 2 < 1^,77) -Rewol < C^) 2 (4.14) 
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for some constant C > independent of y and 77. Therefore if L is large enough and h is 
small enough, we have for v £ Cq°([0, L -1 ]), 

\\{0{t 2 + t + h)(r]) 2 + 0{h)hD t ]v\\ 2 < C- 2 \R{y, 77) - Rew | 2 |M| 2 + 0(/i 2 )||/iAv|| 2 . 

We can apply 

((hD t ) 2 v,v) = \\hD t v\\ 2 + ih 2 D t v(0)^(Q) 

again to get 

\\hD t vf < |((/iD t ) 2 v,«>| + /i 2 |A«(0)||«(0)| < \{\\{hD t fvf + |M| 2 ) + ^1^^(0)^(0)1. 
Therefore 

\\{0{t 2 + t + h)(7]) 2 + 0{h)hD t ]v\\ 2 < C , - 2 (| J R(y, 77) -RewoHMI 2 

+ii(/iA) 2 «ir)+o(^ 4 )iAv(o)iKo)i. (415) 

It is easy to prove the following elementary inequality: 

{Va - h 2 b - VC- 2 a + h%) 2 ^ (1 - 2C- l )a - 2h 2 b (4.16) 
for C large and h small, independent of a, b > 0. 

Applying (4.16) with a = \\\{hD t fvf + {r 2 + \\R{y,r,) - Rew | 2 )|M| 2 , b = (\R(y,rj) - 
w 1 + ^ ro)|Au(0)||i;(0)|, by (4.6), (4. 13) and (4.15), we have 

\\(P(y,t,V, hD t ) - u )v\\ 2 > (1 - 2L7- 1 )[i||(/ i A) 2 ^|| 2 + (r 2 + - FLe 6U 1 2 ) ||^ || 2 ] 

-0(h 2 )(\R(y, V ) - wo | + >/3to)|A«(0)||t;(0)|. 

Now by our assumption, 77) — Rewo| > c = r tan|, and (4.14), we proved the 

lemma. □ 

4.3. Lower bounds for the scaled operator near the boundary. We first consider an 
estimate valid for functions supported in a sufficiently small neighborhood of the boundary. 

Proposition 4.3. Suppose that u £ C oc '(M n \ O) satisfies supp(w) C X x [0, L~ l ), and the 
Robin boundary condition d v u = -fu for 7 £ C°°(X,C). Then 

||(/i 2 P-w H| 2 > \r + S{Reu )hl - 0(h)\ 2 \\u\\ 2 . (4.17) 
Proof. Let T is the FBI transform defined in Section 2.3, by (2.2), we have 
IMIi^Xxfo,!,- 1 )) = II^ m IIl 2 (t*xx[o,l- 1 )) + ^Wll^lli 2 - 

By (2.4), we have 

II M IIh2(Xx[0,L-i)) ~ II( ? ?) 2 ^' U IIl 2 (T*Xx[0,L- 1 )) + II(^A) 2 ^m|| 2 : 2( T , Xx [ 0iL -1)), 
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also 

\\(h 2 P- w )M||ia(xx[o,L-i)) = \\{ p {v^^,hD t ) - ^Tuf^^x^-x^ + 0(h)\\u\\ 
Now (4.2) shows that if \R(y, rj) — Recuol < c, 

POO 

\(P(y,t,rj,hD t ) - uj )Tu\ 2 dt ^ (r + 2S(Reu )h% -0(h)) 2 / \Tu\ 2 dt 



r>oo 

+C- 1 I \(hD t ) 2 Tu\ 2 dt - O(h 2 )\D t Tu(0)\ 2 - O(h 2 )\Tu(0)\ 2 ; 
Jo 

and (4.5) shows that if \R(y, rj) — Reu; | > c, 

poo poo poo 

/ \{P{y,t,r),hD t ) -cu )Tu\ 2 dt^ (r + CT 1 ) 2 / |Tu| 2 di + CT 1 / |(/iA) 2 Tu| 2 rft 
Jo io io 

+C- 1 (r ? ) 4 / |T M | 2 dt-O(/i 2 )(r/) 2 |A^(0)||Tw( 



If we integrate J °° |(-P(y, t,rj, hD t ) — uo)Tu\ 2 dt in (y,rj) e T*X, we get 
II (-P(y, ?7, &A) - w )Tu|| 2 : 2 {r * Xx[ o jr -i )) 

■'l ^(T'lxfO.L- 1 )) 



^ (r + 25(Rew )3/ i 3 - 0(/i)) 2 ||Tu|| 2 



+ C- 1 / / \(hD t ) 2 Tu\ 2 dt+ / (^IT^I 2 ^ 

\Jt*XJ0 J\R(y,rj)-Rcui \>cJo 

- 0{h 2 ) || (r])D t Tu(y, V , 0) ||i 2(r » x) - 0(/* 2 ) || (^>r«(y, 17, 0) || 2 2(T , X) . 
Here D t Tu(y, V ,0) = T(D t u(; 0))(y, 77) = T(-*it(., 0)), so by (2.3), 
|| ( V )D t Tu(y, 77, 0) || 2 2(T » X) =|| (r/)T(fc(-)«(-, 0)) || 2 2(T , X) 

<C7||fc(y)«( yi 0)||^ w < CpH^^H^O)!! 2 ^; 
||(r/)r M (-,0)||l 2(T , x) ^11^,0)11^^. 



(4.18) 



Now we can apply Proposition 2.4 to the last two terms in (4.18) to show that they are 
bounded by 0(/i)||w||#2. 

Notice that if \R(y, rj) — Rewo| < c, then (y, rj) lies in a compact region of T*X, we have 



\R(y,Tj)-Reuj \<c JO 



00 

4|rr„.|2 



<77> 4 |T M |^ ^ C\\Tu\\ 2 L2{T 

•XxfO.L- 1 ))) 



thus 



p poo p poo 

/ / \(hD t ) 2 Tu\ 2 dt+ / / (r7) 4 |T«| 2 d^max{0,C- 1 ||«|| 2 ?2 - C|H|^}. 

JT*XJ0 J\R{y,v)~Rcuj \>c JO 
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Therefore from (4.18), we have 

\\(h 2 P - w )w|| 2 ^(r + 2S(Reu )h% - 0(h)) 2 \\u\\ 2 L2 

+ C7- 1 max{0,C7- 1 || M ||2 2 - C\\u\\ 2 L2 } - 0(h)\\u\\ 2 H2 . 

h h 

This concludes the proof of (4.17). □ 

4.4. Lower bounds for the scaled operator. The estimate away from the boundary is 
now combined with elliptic estimates away from the boundary to give the main estimate 
of the paper: 

Theorem 2. There exists some e > such that for Wo satisfying argwo G (e, | — e), 
Re w G (l-e,l + e), and u G C°°(R n \0) D D(R n \0) (i.e., satisfying the Robin boundary 
condition), we have 

\\(h 2 P - oo )u\\ 2 ^ \r + 2S(Reu )hi - 0(h)\ 2 \\u\\ 2 . (4.19) 

for all sufficiently small h > 0. 



Proof. We only need to estimate the part away from the boundary and connect it with 
(4.17). Let <p ,<pi G C°°(R n ; [0,1]) such that (fil + <ff = 1, suppy? C {x : d(x) < L -1 }, (pi = 
on {x : d(x) > (2L) -1 } where d(x) is the distance from x to O. We claim that 

\\(h 2 P - w )^i«|| 2 > (r + 2S(Recuo)h 2 3) 2 \\(p lU \\ 2 . (4.20) 

In fact, from the argument in Section 2.1, the symbol p of — A|r when d(x) > (2L) _1 takes 
its values in e < — argz < n — e for some e > 0. So by the assumption on w , 

inf \p — wo | > Im(e ie Wo) = |wo| sin(e + argwo) 

2 2 

> | o;o | sin(argwo) + 25*(Re w )3/i 5 
= r + 2S , (Rew )5/ii 

We have the estimate (4.20). 
Since 

{h 2 P — u )(fjU = (fj(h 2 P — w )w — [ipj, h 2 P]u, 

we have 

\\(h 2 P- u )(pju\\ 2 = \\cpj(h 2 P -u )u- [cpj, h 2 P)u\\ 2 

< \\(pj(h 2 P - w )m|| 2 + 2||^(/i 2 P-w )w||||[^,/i 2 P]w|| 
+ \\[<p 3 ,h 2 P]u\\ 2 , 
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thus 

\\(h 2 P-uj )u\\ 2 = J2 \\<Pj(h 2 P -uo)uf 
3=0,1 

^ \\{h 2 p-uj,)^u\\ 2 -Y. fc>* ap Mi a 

- hi(h 2 P-uoM\\[Pi,h 2 P]u\\. 

3=0,1 

The commutators can be estimated by 
|| [^-, /z 2 ^]^!! = <^(/i)(||/i^^^|U 2( ^x [(2 ^ ) -x^- 1D -K II^H) ^ ^(^XH (/i 2 J^ — c^o)^|| + N|) 
since h 2 P — uq is elliptic when d(x) > (2L) _1 . Therefore 

\\(h 2 P - L0 )U\\ 2 > W( h2p ~ W 0)M| 2 - 0(h)(\\(h 2 P - L0 )U\\ 2 + \\ U \\ 2 ). 

3=0,1 

Now we can conclude that 

\\(h 2 P - oo )u\\ 2 0(h)) W( h2p ~ Moduli 2 - 0(h)\\u\\ 2 

3=0,1 

0(h)) \ r o + 2S(Reu )hl - 0(/i)| 2 ||^«|| 2 - 0(h)\\u\\ 2 

3=0,1 

> |r Q + 25 , (Rewb)i/is - 0(/i)| 2 ||w|| 2 . 

□ 

5. The pole-free region 

Now we prove Theorem 1. An equivalent formulation is to say that there are no res- 
onances for P^' in the region Re( > co,0 < — Im£ < S'(ReC) 5 — ci for some constant 
Co, Ci > 0. Suppose £ is a resonance of P (7) such that < -lm( < S(Re()s - c x . Then by 
Proposition 2.1, A = ( 2 is an eigenvalue of P. Let h = (Re^) -1 , then h 2 ( 2 is an eigenvalue 
of h 2 P: h 2 Pu = h 2 ( 2 u for some u G D(R n \ O). Now we apply Theorem 2 to u and 
u = Re(h 2 ( 2 ) + ir ,r > 0. Since hRe( = 1, we have 

Re(h 2 ( 2 ) = h 2 (Re() 2 - h 2 (Im() 2 = 1 + O(hl); 



4 , 



Im(h 2 ( 2 ) = 2/i 2 (ReC)(ImC) = 2Mm( = O(hi). 

So it is easy to choose some r$ such that cuo satisfies the condition in Theorem 2(e.g. r = 1), 
we get 



h 2 C - lu 2 \\\u\\ 2 l2 = \\(h 2 P -u )u\\ 2 L2 ^ \r + 2S(Reu )*h3 - 0(h)\ 2 \\u\ 



L 2 ' 
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where 

h 2 ( 2 -u = i{lm{h 2 ( 2 ) - r ) = i(2hlm( - r„); 
Re w = Re(/i 2 C 2 ) = 1 - h 2 (Im() 2 . 

Thus we have 

|r - 2WmC| 2 ^ |r + 2S[1 - h 2 {Im() 2 ]hl - 0{h)\ 2 . 

so 

-ImC ^ S[l- h 2 {Im() 2 }h-l - M. 

Now by the assumption that — lm( < Sh~a — ci, we can choose cj large, (e.g. ci ^ S' + M), 
so that we have 

1 - h*> (1 - /i 2 (ImC) 2 )5 = (1 - 0(/it))t. 

Let h — > 0, we have a contradiction. Therefore we can choose c , Ci > large such that 
there are no poles in Re( > c , < —ImC < S'(ReC) 5 — C\. 

6. Appendix 

Now we present the proof of Lemma 2.3 following [21]. We need to prove the following 

(1) u extends holomorphically to a function U in a complex open neighborhood of W fl 

(2) ug = U\r g is smooth up to dTg = dO; 

(3) (-A|r fl - X 2 ) ko u e = 0, d a u e \ar e = u a in T g n W. 

Part (1) and the first equation in (3) follows from Lemma 2.2 by choosing intermediate 
contours between Tg and M. n \ O away from the boundary. The difficulty lies in justification 
of the boundary condition for which we need to estimate the norm of Ug. To do this, we 
first review the strong uniqueness property of the scaled operator and its corollary. For 
details, see the appendix of [21]. 

Proposition 6.1. Assume P is an m-th order differential operator with holomorphic coef- 
ficients, T a totally real submanifold of C ra of maximal dimension such that P|r is elliptic. 
Then if u G @'(T) satisfies Pru = on T and u — in a neighborhood of some xo G T, 
then u = 0. 

Corollary 6.2. Let P,T be as in the proposition above, Qi CC ^2 CC ^ C T are open 
sets, then there exists some constant C > such that for all u G H m (Q^), 

\\u\\h^(q 2 ) < Cfll-Pullj^COs) + ||w||flO(n s \ni))- (6-1) 

Also if P, T, fli, 0,2, ^3 depends continuously on some parameters varying in some compact 
set, then we have the estimate for some constant C independent of the parameters. 
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This corollary shows that if Pu = 0, then the part of u in f2 3 \ Qi controls the whole 
part of u in f2 2 - We shall apply this property to a family of intermediate contours between 
Tg and R n \ O and use the part in W 1 \ O to control the part in Tg. 

Now we describe our family of intermediate contours. To do this, we first blow up a 
neighborhood of Xo by introducing the following change of variables: 

x i— > x, x = y + ex 

where y G dO is some boundary point near xq, e > is a parameter which we will let tend 
to 0. We shall choose x such that \x \ = 1 is close to the normal direction to boundary 
through y and focus on the region B(x , 1) in the new coordinates. 

The intermediate contours are constructed as follows: Let Tg ty ^ e be the image of Tg in the 
complexified x-space. Since Tg is parametrized by z = x + i8f'(x), we have the following 
parametrization of ^e,y,e'- 

z = x + iOe' 1 f'(y + ex) = x + i0dxf e , y (x) 

where f e , y {x) = e~ 2 f(y + ex). The derivatives of can be estimated as follows: 

/ 0(|x| 2 ~H) if | a |^2 
OsJcvW | o( e M-2) if \a\>2 ■ 

We choose a cut-off function x £ C£°(B(x , |)),0 ^ x ^ 1 an d X = 1 on -^(^o, !)■ Out 
intermediate contours will be the image of 

x h-> 5 = x + i9d s (xfe, y (x)) 

and we will write ^0,^1,^2,^3 to be the images of the balls B(x , |), B(x , |), -B(xo, |), 
-B(xo, |), respectively. See Figure 2(where we omit fi 2 ). 

By the strong uniqueness property of (— A5 — e 2 A 2 ) fc ° and (6.1), we have the following 
estimates uniformly with respect to y in a neighborhood of xo and e small, 

\M\ H ^o(n 2 ) < C(||(-A 2 - e 2 X 2 ) ko \\ H o m + Hvllflo^no)- 
We shall only use the following weak version: 

IMIz> (ni ) < C{\\{-Ai - e 2 \ 2 ) h °\\ L2m + ||«||^ ( n,\ ni ))- (6-2) 

Since in the x coordinates, u(x) = u(x) satisfies the equation (— A £ — e 2 A 2 ) fc °w = 0. 
By Lemma 2.2, u extends to a holomorphic solution U of (— Aj — e 2 \ 2 ) k °U = over a 
neighborhood of a family of intermediate contours between B(xq, |) and r° ye including 
^2 C ^e,y,e- Back to the x-coordinates, we get the holomorphic extension U of u in W fl r° 
if we let y varies near xo and e goes to 0. 
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Intermediate contours in complex scaling 



Figure 2. Blow up and complex scaling contours. In the complex scaling 
picture, Tg,y,t is the intermediate contour between M. n and ^e,y,e- On this 
contour, Q , Qi, S7 3 are the images of the balls B(x , |), B(x , |), B(x , |), 
respectively. Therefore fl C ^e, y ,e, ^\^iC M. n . 

Substitute v = Un s in (6.2), noticing that Q3 \ Qi C K~, so U = u on Q 3 \ Qi, we have 

II^IU 2 ^!) < C||m||l 2 (c 3 \^i)- 

Back to x-coordinates, we will have similar estimate for U and u with the same constant 
C uniformly for all y near xq and e > 0. In particular, this shows that ug = U\r g is well- 
defined in L 2 near xq. Since ug satisfies the non-characteristic equation (— A z — X 2 )ug = 0, 
if we identify Tg with IR n \ O and use the normal geodesic coordinate (x',x n ) (again, only 
locally near x ), then u e G C([0, e ); ^'(M n_1 )). In particular, it has a boundary value 
u (x', 0) G ^'(W 1 - 1 ). (For the proof of this, see e.g. [9].) 

Now it only remains to show that ug(x',0) coincides with the original boundary value 
u(x'). To do this, we substitute v — U — tt(0) in (6.2), noticing that -u(O) = u(y), 

(-A s - e 2 X 2 ) k °(U - u(0)) = -(-e 2 A 2 ) fco u(0), 
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we have 



U - uWWitW < C((e 2 A) fco |n(0)| Vol(fi 3 )^ + ||« - u(0) Hl^^))- 



Since u is smooth, the last term tends to when e — > 0. Therefore 




e "||w e (x)-u(x / )||| ((:rV)i | ) = o(l),e^0 



(6.3) 



uniformly in x'. Now let Xn £ Co°0^)' x'O^O e C^°(IR n : ) be cut-off functions with support 
close to 1 and 0, respectively. Also let J Xn(x n )dx n = J x'(x')dx' = 1. Then for any test 
function if G C%°(R n - 1 ), 

(u e (x',0) -u(x'),ip(x')) = \im(ue(x) - u(a; / ))e _1 x n (e _1 a; n )^(a; / )da; 



since e n J(ug(x) — u{x'))xn{.£ lx n)x'( e 1 (x' — y'))(p(x')dx has a uniformly compact support 
with respect to y' and tends to zero uniformly in y' by Cauchy-Schwarz inequality and (6.3). 

To get the desired global deformation with suitable boundary condition, we only need to 
glue all the local deformation together using the strong uniqueness property. 

For higher order derivatives, we can repeat the argument for every d a u which satisfies 
the differential equation (—A — X 2 ) k °(d a u) = to get the holomorphic extension of d a u. 
The strong uniqueness property shows that this is exactly the derivative of the holomorphic 
extension of u. 
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